Global properties of the spectrum of the Haldane— Shastry spin chain 
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We derive an exact expression for the partition function of the su(m) Haldane-Shastry spin chain, 
which we use to study the density of levels and the distribution of the spacing between consecutive 
levels. Our computations show that when the number of sites N is large enough the level density is 
Gaussian to a very high degree of approximation. More surprisingly, we also find that the nearest- 
neighbor spacing distribution is not Poissonian, so that this model departs from the typical behavior 
for an integrable system. We show that the cumulative spacing distribution of the model can be 
well approximated by a simple functional law involving only three parameters. 
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I. INTRODUCTION 

The Haldane-Shastry (HS) spin chain describes N 
spins equally spaced on a circle with an interaction in- 
versely proportional to the square of their chord dis- 
tance 1 - 2 . The original motivation for studying this model 
is the fact that it possesses an exact Jastrow-product 
ground state, which coincides with the U — * oo limit 
of Gutzwiller's variational wave function for the Hub- 
bard model^i^, and also with the one-dimensional ver- 
sion of the resonating valence bond state introduced by 
Anderson^. Since its very introduction, the HS spin chain 
has been extensively studied as a completely integrable 
model 7 solvable by the asymptotic Bethe ansatz 8,9,10 , 
whose spinon excitations provide a simple example of a 
system obeying fractional statistics^. 

The energy spectrum of the HS Hamiltonian with spin 
1/2 was partially computed in the original papers of Hal- 
dane and Shastry. In a subsequent publication^, Hal- 
dane et al. empirically found a complete description of 
the spectrum for arbitrary spin, and explained its highly 
degenerate character by the symmetry of the model un- 
der the Yangian algebra y(s\ m ). These results were rig- 
orously established in Ref™ by explicitly constructing a 
transfer matrix in terms of the Dunkl operatorsM^ of 
the trigonometric Sutherland dynamical modeUSii 7 .. In 
this approach, the spectrum is obtained by considering 
all possible motifs S = (OSi . . . 5jy_i0), where each Sj is 
either or 1 and the maximum number of consecutive 
1 's is m — 1 . Indeed, the energy associated with a motif 
5 is given by the compact formula 

JV-l 

E us (S)=^6 j j(j-N). (1) 

The degeneracy of a level Ens is obtained by summing 
the degeneracies corresponding to all the motifs S such 
that Ens(5) = Ens- Although there is a well-defined 
algorithm for computing the degeneracy of each motif, in 
practice the computation becomes quite involved except 
for m = 2. It is therefore difficult to derive in this way 
an exact expression for the partition function valid for 
arbitrary values of N and m. Perhaps as a consequence 



of this fact, little attention has been paid in the literature 
to the global properties of the spectrum of the HS chain. 

Some authors^?' - 9 have suggested that the main obsta- 
cle in computing the partition function of the HS chain 
in closed form is the fact that the dispersion relation ||TJ 
is nonlinear in j, in contrast with the Polychronakos ra- 
tional chainS2i2i. In a recent paper™, however, the parti- 
tion function of the trigonometric HS spin chain of BCn 
type has been exactly computed applying what is known 
as Polychronakos 's freezing trich^, notwithstanding the 
fact that these chains have a nonlinear dispersion rela- 
tion similar to Q). In fact, we shall prove in what follows 
that the partition function of the chain (j2J can also be 
computed using the freezing trick. From the partition 
function it is straightforward to generate the spectrum 
of the HS chain for a wide range of values of N and m, 
and thus study global properties thereof such as the level 
density or the distribution of the spacing between con- 
secutive levels. 



II. PARTITION FUNCTION 

For convenience, we shall take the Hamiltonian of the 
(antifcrromagnetic) Haldane-Shastry spin chain as 

= ij>in(&-e,r 2 (i + %), (2) 

where £j = iir/N and Sij is the spin permutation opera- 
tor of particles i and j. Here and throughout the paper 
all sums and products run from 1 to AT unless otherwise 
specified. The Hamiltonian of the original HS spin chain 
is given by H HS = H - E max , where 

£ max = ^sin(&-^)- 2 (3) 

i<j 

is the highest energy of H. In order to apply the freezing 
trick, we need to introduce the Sutherland spin model 

H * = - E 9 l + a E sin (^ - x ^ 2 ( a + 5 «) > ( 4 ) 

i ijtj 
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and its scalar version 

Hq = — d 2 . + a(a — 1) sin(a;i — Xj)~ 



We thus have 



H* = H + 4aH , 



(5) 



oo 



where H is obtained from H by the replacement £j 
The freezing trick is based on the fact that for a 
the particles "freeze" at the equilibrium positions of the 
scalar part of the potential in H* , which are simply the 
lattice points of the chain @. In this limit, the spin 
degrees of freedom decouple from the dynamical ones, so 
that by Eq. (JSJ the energies of the dynamical spin model 
are approximately given by22 



e: 



O.i 



J > 



(6) 



where Eqa and Ej are any two levels of Hq and H. Hence 
the partition functions Z, Z* , and Zq of H, H* , and Hq, 
respectively, satisfy the approximate equality 

Z*(T) ~ Z (T)Z(g) , a»l. 

The latter equation leads to the exact formula 



Z{T) = lim 



Z*(4aT) 
Z (4aT) 



(7) 



which we will use to compute the partition function of 
the chain (J2J in closed form. 

In order to evaluate the RHS of JJJ, we need to com- 
pute the spectra of H* and of its scalar limit Hq. These 
spectra can be obtained in a unified way by considering 
the scalar differential-difference operator 



H 



i¥=3 



sm{xi-Xj) 2 (a-Py) 



(8) 



where Py permutes the coordinates i and j. The operator 
H is represented by an upper triangular matrix in a (non- 
orthonormal) basis whose elements are of the form 



P (x) 



a 2ipx 



J|sin Q (x l 



(9) 



N is such that 
1, are integers. 



where the vector p = (pi, ■ ■ ■ ,pn) & 
the differences pi — Pi+x, 1 < i < N 
The basis elements (0 should be ordered in a suitable 
way that we shall now describe. We shall say that a 
vector p = (fix, . . . ,Pn) is nonincreasing if p^+i < pi 
for i = 1 , . . . , N — 1. Given two nonincreasing vectors 
p and p', we shall write p -< p' if px — p\ — ■ ■ • = 
Pi—x — Pi-x = and pi < p[. Finally, we say that the 
basis element <f> p precedes <f> p > if p -< p', where p and 
p' are the unique nonincreasing vectors obtained from 
p and p' by reordering their components. It can then 
be shown that the matrix of H in the basis {</> p } with 



the order just defined is indeed upper triangular, with 
diagonal elementsiiSi 



E(p) =J2(2pi + a(N + l-2i)Y 



(10) 



We shall now see how the spectrum of H* follows eas- 
ily from that of H. To this end, let us introduce the total 
antisymmetrizer A with respect to simultaneous permu- 
tations of the spatial and spin coordinates. We can con- 
struct a (non-orthonormal) basis of the Hilbcrt space of 
the Hamiltonian H* with states of the form 



Vw(x) = A(0 p (x)|s)) 



(11) 



where |s) = \sx, ■ ■ ■ , sn) is an element of the spin basis 
and the vector p satisfies the following conditions: 



i) 



Pi+i, 1 < i < N - 1, are 



hi) 



The differences rij = p 
nonnegative integers. 

At most m components of p can be equal. 
The total momentum vanishes, i.e., ^2 i p% = 0. 

The first two conditions are a direct consequence of the 
antisymmetric nature of the states l|llf) . The last con- 
dition reflects the fact that, since H* is translationally 
invariant, we can work in the center of mass frame. The 
basic states il>p, s should be ordered in such a way that 
ip p<s precedes ipp',s' if p -< p' (note that the vectors p 
and p' are nonincreasing by condition i)) 
From the elementary relations PijA = 
fact that H clearly commutes with A, it follows that 

F>p,s=^p,s=A((^ p )|s)) 

= a(p(p)0 p |s)+ Y, c pp'V|s) 
p'^p 

= #(P)'0P,S + E C PP'^P',s ■ 
P'-<P 

Hence the Hamiltonian H* of the Sutherland spin model 
is upper triangular in the basis {^p, s }, with diagonal 
elements 



SijA and the 



E*(p, s)=J2 (2ft + a(N + 1 — 2i)Y 



(12) 



where p satisfies conditions i)-iii) above. 

The spectrum of Hq can be derived by a similar argu- 
ment, noting that Hq = H on scalar symmetric states of 
the form ?p p = A s p , where A s is the symmetrizer with 
respect to the spatial coordinates and p satisfies only con- 
ditions i) and iii) above. Henceii the eigenvalues i?o(p) 
of Hq are also given by the RHS of (|12fl , where now p is 
not restricted by condition ii). 

From the above results it is easy to compute the parti- 
tion functions Zo(4aT) and Z*(AaT) in the limit a — > oo. 
For the computation of Zo(AaT), we start by expanding 
the eigenvalues of Hq in powers of a as 



E (p) = a 2 E° + AaJ2( N + 1 - 2 «)ft + 0(1) 



(13) 



3 



where 



E° = ^2(N + 1 - 2if = i N(N 2 - 1). 



Since E° does not depend on p, and therefore contributes 
the same overall constant factor to both Zo and Z* , we 
shall henceforth drop the first term in Eq. i|13[l . With 
this convention, for d» 1 the denominator in Eq. Q is 
given by 



Z (4ar)~^ g S iPl (JV+i- 



2i) 



where q = e _1 / ( - feBT '' and the outer sum runs over all 
vectors p satisfying conditions i) and iii) above. Setting 
njy = Pn we have 



N-l 



£>(J\T+l-2t) = £>,-(iV+l-2i) = J2 i(N-j)nj . 

i j>i j=l 

Taking into account that tin is determined by the re- 
maining m's by condition iii), we finally obtain 



Z (4aT) 



N-l 

ni,...,njf_i>0 3=1 

N-l 



n(w 

3=1 



j(N-j) 



(14) 



In order to compute the partition function Z*(AaT) for 
a ^ 1, it is convenient to represent the vector p labeling 
the energies (120 of H* as 



(15) 



Note that Yl — N, so that k = (fci, . . . , k r ) belongs to 

i=l 

the set Vn of partitions of N (taking order into account). 
Calling 




Ki — ^ kj , 

3=1 



(16) 



and dropping again the term a 2 E° 7 in the large a limit 
Eq. Q becomes 



A", 



^*(p,s)~4a^ Pl J] (AT+l-2i)=4a^M 

i=l j=Ki-i+l i=l 



where 



l i = k i {N-2K i + k i 



(17) 



Since E*(p, s) does not depend on the spin coordinates s, 
the degeneracy associated with this eigenvalue is given by 



d(k) = n 



so that d(k) = if ki > m for some i, in accordance with 
condition ii). Hence 



Z*{AaT)c d ( k ) E 1- 

kGPjv PX> — >Pr 

fciPlH hfcrPr=0 



(18) 



Calling i/j = p.; — pi + \ £N,i = l,...,r-l, and i/ r = /9 r , 
we have 

r r 

Em= E ^ = E^' ( 19 ) 

1 1^2<j<r / — 1 

where 



N j = Y^l i = K j {N-K j 



i=l 



by Eq. (|17|) . Note, in particular, that the numbers Nj 
depend on k through the partial sums (15)1. Substitut- 
ing l|19|l into l|18[). and taking into account that K r = N 
implies N r = 0, we obtain 



Z*(4aT)^ J2 d ^ E II ^ 
kePiv wi,...,f r _i>0j=i 

r-l 



E n r 



AT, 



(20) 



Combining Eqs. (|14fl and l)20|l. the partition function Z 
can be expressed in closed form as 

AT-1 r-l jy. 

z(T)= n (i E wITr^v-- ( 21 ) 

7=1 ke^jv i=i y 

Note that, by definition, the partial sums Ki are natural 
numbers satisfying 1 < Ki < ••• < i^ r _i < iV — 1. 
Denoting by i"T( <••• < K' N _ r the elements of the set 



and setting 
we have 

JV-l 



N[ = Ki{N-K^) 1 



r-l N — r 

</ v >Hu <r 



n (i-^)) =n( i-^)n(i v 

j = l i=l i=l 



This identity and Eq. i|21|) yield the following remarkable 
formula for the partition function of the spin chain @ : 

r / \ r— * N — r 
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From the previous formula it follows that the energy 
levels of H are of the form 



N-l 



(23) 



where 5j = 1 if j is one of the partial sums Ki corre- 
sponding to a partition (fci, . . . ,k r ) £ Vn with fc; < m 
for all I (by condition ii)), and Sj = otherwise. In or- 
der to relate Eq. lll'-'it with the known expression (|TJ for 
the energies of the original HS Hamiltonian, we need to 
evaluate the maximum energy E max . From Eq. (Q we 
have 



EL 



N-l 

2 f ( .7 )^ 



,2 ^ 

.TV 



= E 



J CSC 



I>» 2 (f 



Hence 



iV-l 



^ = fE- 2 (f) = ^ 2 -i) 



(24) 



3=1 



where the last sum is evaluated in RefiS. Since the RHS 
of coincides with the sum J2f=i j( N - j)i Ec l- I23> 
implies Eq. 0) with <5j = 1 — In particular, from the 
latter relation between S and 5' it follows that S is a motif 
with no more than m — 1 consecutive l's. 



III. LEVEL DENSITY AND SPACINGS 
DISTRIBUTION 

The RHS of Eq. (|22|l is a polynomial in g whose evalu- 
ation with a symbolic algebra package is straightforward 
once N and m are fixed. In this way we have been able to 
compute the spectrum of the chain @ for relatively large 
values of N and m, for which the usual motif approach 
becomes inefficient due to the difficulty of computing the 
degeneracies. From the analysis of the spectral data thus 
obtained one can infer several global properties of the 
spectrum that we shall now discuss. In the first place, it 
is apparent that for N 3> 1 the level density is Gaussian 
to a very high degree of accuracy, as in the HS spin chain 
of BCat type studied in Refi^S. In other words, for large 
iV the cumulative level density 



F{E) 



is approximately given by 



-A' 



E < 

i;Ei<E 



G{E) = \ 



1 + erf 



E-fi 

V2<7 



where is the degeneracy of the energy Ei, and \i and 
cr are respectively the mean and the standard deviation 
of the energy. This can already be seen, for instance, 
in the case N = 15 and m = 2 presented in Fig. ^ 
The agreement between F and G rapidly improves as N 
and/or m grow, e.g., for m = 2 the mean square error 
decreases from 5.2 x 10~ 5 for N = 15 to 5.6 x 10 -6 for 
N = 20, or from 2.6 x 1(T 5 for TV = 15 to 2.6 x 1(T 6 for 
N ~ 20 when m = 3. 

F(E), G(E) 




FIG. 1: Cumulative distribution functions F{E) (at its dis- 
continuity points) and G(E) (continuous line) for iV = 15 and 
m = 2. 

Since, by the previous discussion, for large N the level 
density is characterized by /i and a through the Gaussian 
law, it is of interest to compute these parameters in closed 
form as functions of N and m. In the first place, using 



the identity trSV,- 
obtain 



m N - 1 and Eqs. © and pijl. we 



1 



2m 



csc 2 (& - £j 



m + 1 
12m 



N(N 2 - 1). 



Similarly, the formula 



tr(SijSki 



= m N-2+2S tk S jl +28 il S jk 



yields 



tr(H 2 ) (tiH) 



,2N 



Am 2 



csc 4 (& - £/) 



(m/ 



-1)N N ^ 



8m 2 

m 2 — 1 
360m 2 



esc £j 
N(N 2 - 1)(N 2 



11) 



(cf. Refi^i for the last equality) . 

The level density is also Gaussian as iV — > oo for the 
so-called "embedded Gaussian ensemble" (EGOE)2& in 
Random Matrix Theory. Note, however, that in the 
EGOE this property is valid provided that the number of 
one-particle states tends to infinity faster than N. This 
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additional condition clearly does not hold in our case, 
since the number of one-particle states (i.e., to) is fixed. 
Another characteristic feature of the EGOE is the fact 
that the nearest-neighbor spacing distribution p(s) is ap- 
proximately given by Wigner's law 

p(s) = (7r/2)sexp(-7rs 2 /4), 

as for the classical Gaussian orthogonal ensemble 27 . On 
the other hand, since the HS spin chain is integrable, 
one would expect that its nearest-neighbor spacing dis- 
tribution obey Poisson's lawp(s) = c~ s , according to the 
conjecture of Berry and Tabor for a generic integrable 
model 2 ^. This conjecture has been verified for a variety 
of integrable many-body problems, such as the Heisen- 
berg chain, the t-J model, the Hubbard model 2 ^, and 
the chiral Potts models. One of the main results of this 
paper is the fact that the nearest- neighbor spacing distri- 
bution of the HS chain deviates substantially from both 
Wigner's and Poisson's laws. 

In order to correctly take into account the effect of 
the local level density in the study of p(s), one must 
first apply to the "raw" spectrum the so-called unfold- 
ing mapping^!. This mapping is defined by decomposing 
the cumulative level density F(E) as the sum of a fluc- 
tuating part Ffi(E) and a continuous part which 
is then used to transform each energy Ei, i — 1, ... ,71, 
into an unfolded energy = £(Ei). The function p(s) 
is defined as the density of the normalized spacings 
Si = (6+1 - where A = (£„ - - 1) is the 

mean spacing of the unfolded energies. By the previous 
discussion, in our case we can take the unfolding map- 
ping £ (E) as the cumulative Gaussian distribution G(E) 
with parameters fi and a given by the previous formu- 
las. As for the level density, to compare the discrete 
distribution function p(s) with a continuous distribution 
it is more convenient to work with the cumulative spac- 
ing distribution P(s) — J Q p(x)dx. Our computations 
for a wide range of values of N and m show that P(s) 
is essentially different from either Poisson's or Wigner's 
law, since its slope tends to infinity both as s — > and 
s —> Smax, where s max is the largest spacing. In fact, it 
turns out that in all cases P(s) is well approximated by 
a cumulative distribution of the simple form 



P(s) 



p(s) = t a [i- 7 (i-ty 



(25) 



where t = s/s max and < a, /3 < 1. The parameter 7 is 
fixed by requiring that the average spacing be equal to I , 
with the result 
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f 



f 



)/B(a+l,P+l), (26) 



where B is Euler's Beta function. For instance, for N = 
26 and m = 2 the largest spacing is s max = 3.06, and the 
best least-squares fit parameters a and j3 are respectively 
0.31 and 0.23, with a mean square error of 4.1 x 10~ 4 (see 
Fig. EJ. 




2.5 3 



FIG. 2: Cumulative spacing distribution P(s) and its approx- 
imation P(s) (grey line) for N = 26 and m = 2. For conve- 
nience, we have also represented Poisson's (long dashes) and 
Wigner's (short dashes) cumulative distributions. 



For a fixed value of m, the parameters a, (3 and s max 
vary smoothly with N > 15, provided that N has a fixed 
parity^. For instance, in Fig.|3|we plot these parameters 
for Tn — 2 and odd N running from 15 to 27 (the plot for 
even N is very similar). In all cases, the fit of the distri- 
bution (|25|l to the data is quite good, the mean square 
error never exceeding 7.4 x 10 -4 . We have performed a 
similar analysis for m = 3 and 15 < N < 22, obtaining 
totally analogous results. 



A' 



FIG. 3: Values of a (box), /3 (rhombus), and s n 
for m = 2 and odd N. 



t /10 (cross) 



The divergence of the nearest-neighbor spacing distri- 
bution p(s) for small s is probably related to the flatness 
of the tail of the Gaussian distribution. It could also 
be argued that, since the Haldane-Shastry chain is com- 
pletely integrable, the full spectrum is a superposition of 
the spectra of the Hamiltonian restricted to subspaces of 
common eigenfunctions of a suitable family of commut- 
ing first integrals. It is well known, in this respect, that 
a superposition of a large number of unrelated spectra 
leads to a sharp increase in the number of very small 
spacings 32 . On the other hand, we do not have a clear 
explanation of the fact that p(s) also diverges when s 
approaches the largest spacing s max . This fact, which 
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P(s) 



0.25 0.5 0.75 1 1.25 1.5 1.75 

FIG. 4: Cumulative spacing distribution P(s) (solid dots) for 
states with zero total spin and odd parity when N = 13 and 
m — 2. For comparison purposes, we have represented Pois- 
son's (long dashes) and Wigner's (short dashes) cumulative 
distributions. 



certainly deserves further study, could be a characteris- 
tic property of all spin chains of Haldane-Shastry type. 

Our results also imply that Berry and Tabor's con- 
jecture does not hold for the HS spin chain, even if we 
restrict ourselves to a subspace of the whole Hilbert space 
with well-defined quantum numbers. Indeed, the nearest- 
neighbor spacing distribution of the superposition of even 
a small number of spectra with Poisson-distributed spac- 
ings must also be of Poisson typa^. As an illustration of 
these assertions, we present in Fig.0]a plot of the cumula- 



tive spacing distribution corresponding to the restriction 
of the Hamiltonian J2J to the subspace with zero total 
spin and odd parity for N = 13 and m = 2, obtained by 
a numerical computation of the spectrum of H restricted 
to this subspace. It is apparent from this plot that P(s) 
is neither Poissonian nor of Wigner type, and that it 
is well approximated by a function of the form i|25|) for 
spacings s > 0.25. It is also clearly noticeable that p(s) 
tends to infinity as s approaches the maximum spacing 

^max — 1-73. 

The non-Poissonian behavior of the spacing distribu- 
tion could in principle be due to finite-size effects^ 3 .. Al- 
though this possibility should be explored in more detail, 
our data clearly show that the cumulative spacing distri- 
bution P(s) is of the form (|25|l for a wide range of values 
of N < 27. 

Note, finally, that an interesting integrable model not 
obeying the Berry-Tabor conjecture has been recently 
constructed in Ref. 34 . In contrast with the HS spin chain, 
the latter model is a non-generic element of a class de- 
pending on a large number of parameters, and involves 
many-body interactions. 
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